Modelling ignition probability with pairwise mixing-reaction model for flame particle tracking

Reduced order models for ignition analysis can offer insights into ignition processes and facilitate the combustor optimization. In this study, a Pairwise Mixing-Reaction (PMR) model is formulated to model the interaction between the flame particle and the surrounding cell mixture during Lagrangian flame particle tracking. Specifically, the model accounts for the two-way coupling of mass and energy between the flame particle and the surrounding shell layer by modelling the corresponding turbulent mixing, chemical reaction and evaporation process if present. The state of a flame particle, e.g., burnt, hot gas or extinguished, is determined based on particle temperature. This model can properly describe the ignition process with a spark kernel being initiated in a nonflammable region, which is of practical importance in certain turbine engines and has not been rigorously accounted for by the existing models based on the estimation of local Karlovitz number. The model is integrated into an ignition probability analysis platform and is demonstrated for a methane/air bluff-body flame with the flow and fuel/air mixing characteristics being extracted from a non-reacting simulation. The results show that for the spark location being at the extreme fuel-lean outer shear layer of the recirculation zone, PMR can yield ignition events with a significant number of active flame particles. The mechanisms for the survival of the initial flame particles and the entrainment of the survived flame particles into the recirculation zone are analyzed. The results also show that the ignition probability map from PMR agrees well with the experimental observation: a high ignition probability in the shear layer of the recirculation zone near the mean stoichiometric surface, and low ignition probabilities inside the recirculation zone and the top stagnation region of the recirculation zone. The parametric study shows that the predicted shape of the ignition progress factor and ignition probability is in general insensitive to the model parameters and the model is adequate for quantifying the regions with high ignition probabilities.     

基于DQN的反舰导弹火力分配方法研究

反舰导弹作为海上作战的主战武器，由于其精度高、射程远、威力大等特性长期以来一直被当作舰艇编队的主要防御对象。针对反舰导弹打击舰艇编队的火力分配问题，我们提出了一种基于深度Q值网络求解反舰导弹火力分配策略的算法。不同于现有的基于领域知识的方法，深度Q值网络无需依赖任何先验信息，就能够通过与模拟器的交互自动求解最佳的攻击策略。该算法使用深度神经网络拟合Q值函数，解决了传统强化学习中的状态空间过大无法遍历的问题。实验结果表明，在各种不同的舰队防御配置下，深度Q值网络求解得到的攻击策略均获得了最佳的毁伤效果。     

Orbital stability close to asteroid 624 Hektor using the polyhedral model

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid’s mass center when the density increases, and the inner equilibrium point moves close to the asteroid’s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary’s equatorial plane of this binary Trojan asteroid; the results indicate that the orbit remains stable after 28.8375?d.     

全模颤振风洞试验三索悬挂系统多体动力学分析

全模颤振风洞试验需要通过软支撑系统模拟飞行器的自由飞行状态并调整模型姿态达到配平状态。参考NASA双索悬挂方案，提出了一种两电机驱动的三索悬挂系统，利用后方两索的同向/反向联动实现模型俯仰和滚转姿态的调整，利用弹簧刚度以及钢绳张力设计实现支撑频率要求。基于柔性多体动力学方法，建立了包括飞行器刚体模型、柔性索、滑轮、弹簧、气动力模型、伺服电机控制在内的复杂系统动力学模型，其中，利用任意拉格朗日-欧拉（ALE）变长度索单元描述钢绳，利用不约束物质坐标的索结点约束描述钢绳与滑轮相互作用，利用索结点物质输运速度约束描述伺服电机绞盘，利用飞行力学的气动力模型描述吹风下的气动力。基于该模型，通过小扰动响应辨识研究了弹簧刚度、钢绳张力、连接点位置等因素对支撑频率的影响规律，并分析了系统姿态调整能力，俯仰调整范围达到-12.5°~12.5°，滚转调整范围达到-45°~45°。采用滑轮处电位计测量的钢绳相对位移作为反馈信号，基于设计的控制律利用多体动力学求解器与Simulink对风洞吹风下的姿态调整过程进行仿真，模型达到配平状态，获得了吹风下的索拉力和伺服电机功率，为系统设计提供基础。     

电动式磁悬浮动量球电磁场耦合分析

本文针对电动式磁悬浮动量球(Electrodynamic suspension reaction sphere, EDSRS)的电磁场耦合进行分析，采用定子绕组的感应电动势评判耦合强度，提出了耦合强度评判指标，作为EDSRS的结构设计约束，以降低电磁场耦合的影响。首先，结合EDSRS的具体结构特点，对多定子之间的电磁场耦合问题进行分析，采用定子电磁场在另一个定子绕组中产生的感应电动势评判耦合强度。然后，对感应电动势产生规律进行分析，提出了耦合强度评判指标，并基于所提评判指标，建立了耦合约束验证流程，作为EDSRS的结构设计约束，最后，通过有限元仿真和实验验证了上述电磁场耦合分析。     

基于深度学习的空间站舱内服务机器人视觉跟踪

为提升舱内跟随服务机器人的任务辅助能力,解决机器人对航天员的视觉跟踪问题,提出了一种基于深度学习和概率模型的人体视觉跟踪算法。利用深度卷积神经网络实现了对穿着多样、姿态任意人体的稳定检测。结合人体检测结果,设计了运动预测概率模型,实现了对指定人员准确、连续的跟踪。算法对包含大多数航天员活动的多个数据集进行了验证。实验结果表明:提出的跟踪算法实现了对穿着多样、姿态任意人体的稳定跟踪,并有效避免了由于穿着相似、遮挡可能造成的误跟踪问题。该算法为空间站舱内跟随服务机器人对航天员的视觉跟踪提供了有效的解决方法。算法基于融合的RGB-D图像,工程上易于构建和实现,也可拓展到其他跟随服务机器人视觉跟踪任务中。     

新一代战斗机非定常流动数值研究综述

新一代战斗机强调超机动能力和强隐身性，其中大攻角下的静态失速、动态失速及内埋弹仓绕流是与高机动和强隐身密切相关的、极具挑战的几类典型的非定常流动，它们对数值方法提出了极高的要求。为了高精度地仿真流场、清楚地揭示流动机理，有效地控制非定常流动，非常有必要发展高精度且高效率的RANS-LES混合方法体系，包含RANS-LES混合方法本身、与RANS-LES混合方法匹配的高精度自适应耗散格式、基准湍流模式、高质量计算网格、高精度时间推进方法、非定常量的统计方法等，具有极强的紧迫性。提出、发展、验证并应用该类方法数值仿真新一代战斗机（包括单独部件、组合体、甚至全机）的非定常流动，数值预测结果与风洞实验数据吻合良好；此类方法可为新型战斗机设计提供理论依据和分析手段。     

空间碎片能量转化利用技术

针对目前空间碎片问题，提出空间碎片发动机概念，立足于使用捕获到的空间碎片，转化为发动机可用的推进剂。在完成碎片清理目标的同时，获得可持续的动力来源，延长清理器的工作寿命。针对空间碎片制粉的方法进行研究，提出使用球磨仪对金属样本进行研磨。使用转刀式粉碎机对非金属材料进行粉碎。通过实验发现，多数粉末粒径达到微米量级。针对空间碎片粉末推进方式进行研究，提出使用静电加速推进方式对粉末进行加速。空间碎片发动机虽然起源于空间碎片清理任务，但是可持续的推进剂供应，也将为小行星探测等任务提供更好的思路。     

Generalized Kutta–Joukowski theorem for multi-vortex and multi-airfoil flow with vortex production——A general model

By using a special momentum approach and with the help of interchange between singularity velocity and induced flow velocity, we derive in a physical way explicit force formulas for twodimensional inviscid flow involving multiple bound and free vortices, multiple airfoils, and vortex production. These force formulas hold individually for each airfoil thus allowing for force decomposition, and the contributions to forces from singularities(such as bound and image vortices,sources, and doublets) and bodies out of an airfoil are related to their induced velocities at the locations of singularities inside this airfoil. The force contribution due to vortex production is related to the vortex production rate and the distance between each pair of vortices in production, thus frameindependent. The formulas are validated against a number of standard problems. These force formulas, which generalize the classic Kutta–Joukowski theorem(for a single bound vortex) and the recent generalized Lagally theorem(for problems without a bound vortex and vortex production) to more general cases, can be used to identify or understand the roles of outside vortices and bodies on the forces of the actual body, optimize arrangement of outside vortices and bodies for force enhancement or reduction, and derive analytical force formulas once the flow field is given or known.     

Drag Reducing and Increasing Mechanism on Triangular Riblet Surface

Drag reducing and increasing mechanism on riblet surface has been studied through computational fluid dynamics(CFD).Drag reduction is achieved through the optimization of riblet geometry which would affect flow structure inside riblet grooves.Force and flow structure on riblet surface are analyzed and compared with those of smooth surface based on the k-εturbulence model.Drag reducing and increasing mechanism is proved to be related to microvortexes induced inside riblets which lead to Reynolds shear stress reduction significantly and is considered to be the dominant factor resulting in wall friction reduction.Simulation results also show that the pressure drag generating from the deviation of static pressure on the front and rear ends of riblets occurs and grows exponentially with Mach number,which can cause drag increasing.Furthermore,near-wall vortical structures,Reynolds shear stress and static pressure on riblet surfaces are also analyzed in detail.